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The method of generating a family of new solutions starting from any wave function satisfying 
the nonlinear Schrodinger equation in a harmonic potential proposed recently [ J. J. Garci'a-Ripoll, 
V. M. Perez-Garci'a, and V. Vekslerchik, Phys. Rev. E 64, 056602 (2001)] is extended to many- 
body theory of mutually interacting particles. Our method is based on a generalization of the 
displacement operator known in quantum optics and results in the separation of the center of mass 
motion from the internal dynamics of many-body systems. The center of mass motion is analyzed 
for an anisotropic rotating trap and a region of instability for intermediate rotational velocities is 
predicted. 
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I. INTRODUCTION 



In a recent paper Garcfa-Ripoll, Perez-Garci'a, and 
Vekslerchik Q have studied the properties of the solu- 
tion of the nonlinear Schrodinger equation in a time- 
dependent, anisotropic harmonic potential 

id t i>(r,t) = (~A + \r-A{t)-v + G(\^(r,t)\)\ </>(r,i), 



(1) 



where A(t) is an arbitrary symmetric positive 3x3 ma- 
trix and G is an arbitrary real function. They have shown 
that from every solution of this equation one may obtain 
the whole family of new solutions by a translation ac- 
companied by a change of the phase 



ip(r, t) V'(r, t) = tl>(r - R(t), t)e ie(r ^\ 



(2) 



provided R(i) is a solution of the classical equation of 
motion in the harmonic potential 



d 2 R(t) 
dt 2 



= -A(t>R(t) 



and the phase 6(r, t) is given by the formula 



(3) 



(4) 



where f(t) is the classical action calculated along the 
trajectory R(t) (a factor of —1/2 is missing in Eq. (12) 
of Ref. fi) 



1 fdR(t) dH(t) 



/m = 3 M 



dt 



dt 



B.(t)-A(t)-K(t) . (5) 



This result implies that the motion of the center of the 
wave packet is governed by the classical equation (^) and 
it decouples from the internal motion that determines 
the shape of the wave packet. It has been argued in 
that such a decoupling may be of significance for the 
dynamics of Bose-Einstein condensates, since they are 
often described in the mean-field approximation by the 
nonlinear Schrodinger equation. 

In this paper we extend the conclusions of Ref. Q| 
to the general case of a mutually interacting many-body 
quantum system. Our generalization is inspired by the 
observation that the new solutions (|^) may be generated 
from the initial solution by the displacement operator 
well known in quantum optics p, vm 



i(r-P(t)-p-R(t)) 



(6) 



where f and p = — zV are the quantum mechanical po- 
sition and momentum operators whereas R(f) and P(t) 
represent the phase-space trajectory of a classical parti- 
cle in the harmonic potential. In order to show that the 
transformation 



^(r,i)^V'(M) 



^(r-P(t)-p-R(t)) 



iP(v,t) 



(7) 
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is the same as that given by (Q), we apply the Baker- 
Hausdorff formula (using the canonical commutation re- 



2 



lations between f and p) to write the displacement oper- 
ator in the form 



»(f.P(t)-p-R(t)) 



-iR(t)-P(t)/2-ir-P(t)_-ip-R(i) 



(8) 



The last term on the r.h.s shifts the argument of the wave 
function by — R(£) and the middle term produces the first 
term of the phase @). Finally, we note that the action 
f{t) after integration by parts reduces to the boundary 
terms 



/(^«.R(t) 

3 K ' 2 dt w 



(9) 



Therefore, the first term on the r.h.s. of (||) reproduces 
the second term of the phase ([|) (up to an irrelevant 
constant phase given by the lower boundary term). 

The action of the displacement operator on a solution 
of a nonlinear Schrodinger equation leads to a shift of the 
center of the wave packet. The vector R(i) that deter- 
mines this shift is a solution of the classical equations of 
motion in the harmonic potential. Therefore, the center 
of the wave packet always follows a classical trajectory. 
In the present paper, we cary out a similar analysis for an 
interacting many-body system. The role of the center of 
the wave packet is played by the center of mass operator 
of the whole system. 

In Section II we construct the transformation analo- 
gous to (0) in the many-body theory of a system of mu- 
tually interacting bosons (or fermions) and we show that 
it also transforms solutions of the Schrodinger equation 
into new solutions of this equation. In Section III we 
show that all these properties are a result of a complete 
decoupling of the center of mass dynamics from the in- 
ternal dynamics of the system confined in a harmonic 
potential. This decoupling has been noticed before (cf., 
for example, Ref. Q). In Section IV we analyze the spec- 
trum of stationary states for a rotating trap and we point 
out that the decoupling causes the splitting of the spec- 
trum. We also show that for an anisotropic trap there 
always exists a region of instability at intermediate rota- 
tional velocities. 



II. DISPLACEMENT OPERATOR IN 
MANY-BODY THEORY 



The many-body theory of interacting atoms will be 
described within the formalism of second quantization. 
We shall use systematically the Schrodinger picture - 
all operators will be time independent. The Hamiltonian 
of the system of atoms contained in a harmonic trap (in 
natural units, m = 1, h = 1) has the form 

H(t) = i fd 3 rft(r) (-V- V + r-i(f)-r) ip(r) 
+ i fd 3 r fd 3 r'ft(r)ft(r r )V{r-r')ifj{r r )iP(r), (10) 



where A(t) is the matrix of the harmonic potential of the 
previous Section and V (r— r') is an arbitrary two-particle 
interaction potential. The field operators ifi(r) and ^(r) 
that annihilate and create particles at the point r obey 
the standard Bose-Einstein or Fermi-Dirac commutation 
relations (we disregard the spin degrees of freedom) 



^(r),^!-')! =<5 (3) (r-r'). 



(11) 



In the construction of the displacement operator in many- 
body theory (following our earlier works |^] ) we employ 
the operators of the total number of particles N, of the 
of the total position R, and of the total momentum P: 



N = d 3 rft(r)4>(r), 



(12a) 

R = fd 3 rft (r)r^(r), (12b) 
P = -i [d 3 r$(r)V4>(r). (12c) 



The number of particles operator N is a constant of mo- 
tion and it commutes with R and P. The operators R 
and P satisfy the commutation relations 



R k ,Pi 



iNSm. 



(13) 



We define the many-body unitary displacement opera- 
tor D(t) by the same general formula (|J) but with second- 
quantized operators R and P playing now the role of 
position and momentum, 



D(t) 



_ „i(R-P(t)-P-R(t)) 



(14) 



We would like to note that the new displacement opera- 
tor D(t) does not share all the properties of the Glauber 
displacement operator in quantum optics. Namely, it 
does not generate a coherent state when acting on the 
vacuum state but leaves the vacuum state unchanged 
D(t)\0) = |0). This is due to the fact that the many- 
body vacuum state is not an exact counterpart of the 
QED vacuum. 

In the proof of the invariance of the Schrodinger equa- 
tion 



WtMt)) = H(t)Mt)) 



(15) 



under the action of D(t) we shall again employ the Baker- 
Hausdorff decomposition of the displacement operator: 

e i(RP(t)-p.R(t)) = e - J JVC(t) e »RP(t) e -»P-R(t) j ( 16 ) 

where C{t) = R(t)-P(t)/2. Next we observe that the 
three factors on the r.h.s. of this formula act in the fol- 
lowing way on the annihilation operators 



e iNC{t)^ r y-iNC{t) = e -iC(t)^ (r)) 
e -iR.P(t)^ r)e iR.P(t) = e ir-P(t)^( r)) 
e iP.R(t)^ r ) e -iP.R(t) = ^ r _ R ^) ) 



(17a) 
(17b) 
(17c) 
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where we have used the following commutation relations 
between the operators N, R, and P and the annihilation 
operators, 



AT,V(r)J = -V(r), 
R.^r) = -rV>(r), 
P,V>(r) = iVi>(r). 



(18a) 
(18b) 
(18c) 



The formulas for the creation operators are obtained by 
hermitian conjugation. With the help of Eqs. ( |l7| ) we 
may calculate the action of the displacement operator on 
the Hamiltonian 



D\t)H{t)D{t) = -jd 3 rft(r) (-iV + P(*)) t/>(r) 
+ i Jd 3 r^(r){r + R{t))-A(t)-(r + R(t))ifj(r) + ^ Jd 3 r Jd 3 r'ft (r)^ (r')V(r - r')^(r')V>(r) 
= + P-P(t) + R-A(t)-K{t) + — (P(t)-P(t) + R{t)-A(t)-R(t) 



(19) 



r 



The interaction term in the Hamiltonian does not change 
because it is invariant under translation and under the 
change of the phase of the field operators. The last ingre- 
dient needed to prove the invariance of the Schrodinger 
equation is the following transformation formula for the 
time derivative 



b\t)id t b{t) = id t 



N d(R(t)-P(t)) 



dt 



R 



dP(t) - dR(t) 



dt 



dt 



NTL(t) 



dP{t) 
dt 



(20) 



All the terms in this formula that contain functions R(<) 
and P(t) cancel with their counterparts in the trans- 
formed Hamiltonian (with the use of the classical equa- 
tions of motion) and we finally obtain 



id t D(t)Mt))=H(t)D(t)\*(t)). 



(21) 



Therefore, for each classical trajectory we obtain a new 
solution of the many-body Schrodinger equation, as was 
the case for the nonlinear Schrodinger equation. 

In order to connect this result with the effective one- 
particle theory described in the Introduction, we shall 
represent the state vectors by the n-particle wave func- 
tions ^(ri,r 2 , . . . ,r„,i), 



*(n, r 2 , . . . , r n , i) = (0|V(ri)V(r 2 ) • ■ • *A(r„)|*(i)>. (22) 



Using the unitarity of the displacement operator and the 
invariance of vacuum state, we obtain (with the help of 
the formulas (|l7|)) the following relation between the the 
wave functions of the initial ^(t)) and the transformed 
state D(t)\V(t)) 



*(n,r 2 ,...,r„,t) -» *'(n,r 2 ,...,r„,i) = (0|^(n)^(r 2 ) . . . j>(r n )D(t)\* (t)) 

= (o|^(t)^t(t)^( ri )^(t)^t( t )^( ra )^( t )^t( t ) ...... b{t)b\t)4>{v n )b{t)\y{t)) 

= e <9(n 1 t)-H9(r 2 ,t)+...-Hfl(r n ,t)^( ri _ R( t ), ra - R(t), . . . ,r„ - R(t), t), (23) 



where the phase 0(r,t) is given by the same formula (Q) 
as in the one-particle theory. Thus, to each solution of 
the many-body Schrodinger equation there corresponds 
the whole family of solutions labeled by the classical tra- 
jectories in a harmonic trap. The new many-body wave 
functions are obtained from the original wave function 
by exactly the same transformation as in the one-particle 



case: a time-dependent shift of its arguments and a multi- 
plication by the same phase factors. The transformation 
formula (E3) acquires a particularly simple form in the 
special case, when the initial wave function describes a 
state in which the center of mass is not entangled with 
the internal degrees of freedom, 

#(ri, r 2j . . . , r„, t) = *(p, £ 2 , ■ ■ • , £n-i, t), (24) 
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where p = (r x + r 2 H h r n )/n and £i,£ 2 , ■ • -,£n-i 

are some relative coordinates, invariant under the trans- 
lations Vi — > Vi + a. Then, the transformation (^) affects 
only the center of mass wave function 



*(P,t) 



$(p-R(t),t), 



(25) 



leaving the wave function of the internal motion intact. A 
general many-body wave function can always be written 
as a sum of product wave functions (|24|). The action 
of the displacement operator on a general wave function 
describing a state in which the center of mass motion is 
entangled with the internal motion gives 



£(t)*(ri,r 2 ,...,r n ,t) = e l8 ^ 
X £fc*k(p- R-(*),*)**(€i> €2. • ••.€»- 



,t) . (26) 



In this case the integration over the center of mass posi- 
tion (either p or R(t)) introduced in Ref. 0] leads to a 
mixed internal state of the condensate. The results pre- 
sented in this Section are closely related to the fact that, 
as is shown in the next Section, the total Hamiltonian 
( |l0| ) may be split into the Hamiltonian of the center of 
mass and the Hamiltonian of the internal motion. 



III. DECOUPLING OF THE CENTER OF MASS 
MOTION 

The separability of the center of mass dynamics from 
the internal dynamics in a harmonic trap follows directly 
from the form ( |To|) of the Hamiltonian. Indeed, this 
Hamiltonian may be written as a sum of two commut- 
ing parts: the Hamiltonian of the center of mass motion 
Hcwi(~t) and the Hamiltonian H\{t) describing the inter- 
nal structure of the condensate, 



H{t) = ff C M(t) + i?iW, 
PP R-A(t)-R 



Hqm (t) 



2N 



2N 



(27) 
(28) 



Hi(t) = \ J<?rft{r)(- VcrVc + r c -A{t)- r c )$(r) 
+- Jd 3 rjd 3 r'^ (r)^(r')y(r - r')^(r')^(r), (29) 

where Vc and rc are the one-particle operators shifted 
by the (normalized) center of mass operators 



r c =r-R/iV, -iVo = -iV-P/iV. 



(30) 



The sum ( |27| ) describes the dynamics of two independent 
systems since the center of mass operators R and P com- 
mute with the Hamiltonian Hi (t) , despite the appearance 
of these operators in H\(t). To prove this statement we 
observe that 



R, Hi(t) = R, H(t) ~ R,H C M(t) =P-P = 0. (31 



Similarly one can show that P, Hj(t) =0. 

The many-body displacement operator (|l4|), con- 
structed from the center of mass operators R and P, com- 
mutes with Hi(t). Therefore, the transformation gener- 
ated by the displacement operator does not change the 
internal state of the system; it only acts on the center of 
mass variables. 



IV. SEPARABILITY OF THE SPECTRUM 



The Hamiltonian ( pL 0| ) is time dependent because we 
have allowed for a time dependent trap potential. Such 
a Hamiltonian does not possess true stationary states. 
The most interesting case of a time-dependent harmonic 
potential occurs when the time dependence is caused by 
a rotation of the trap. Then, by going to the rotating 
frame we can eliminate the time dependence and study 
the spectral properties of the resulting Hamiltonian. The 
price to be paid for this simplification is the appearance 
of an extra term in the Hamiltonian — the scalar product 

r2 M of the angular velocity H and the total angular 
momentum of the system M, 



M = -i /d 3 r^ + (r)(r x V)^(r). 



(32) 



This term is responsible for the inertial forces: the Corio- 
lis force and the centrifugal force. The total Hamiltonian 
in the rotating frame is 



Hr = HcM+H h 
- P P R -A R RxP 

ACM = 1 *l ' - 



2N 



2N 



N 



(33) 
(34) 



#1 = 5 fd 3 rft(r) (-Vc-Vc +r c -i-r c ) ^(r) 
+ - fd 3 rfd 3 r'^(r)^(r')V(r~ r')^(r')^(r) 



+ ifl- d 3 rft (r)(r c x V c )^(r), 



(35) 



where A is the time independent matrix whose eigen- 
values a x , a y , and a z are the squared frequencies of the 
trap potential. We use the coordinate system with the 
axes directed along the principal directions of the trap 
potential. 

The transformation to the rotating frame preserves the 
separability of the center of mass motion from the inter- 
nal motion. Therefore, the spectrum of the many-body 
Hamiltonian is a Cartesian product of the center of mass 
spectrum and the spectrum of the internal motion. Each 
eigenvalue of the Hamiltonian Hi gives rise to the whole 
ladder made of the levels of the center of mass Hamil- 
tonian Hcm- These levels can be exactly calculated re- 
gardless of the form of the mutual interaction. Since for a 
fixed number of particles the commutation relations be- 
tween R and P are canonical (up to a numerical factor 
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of n), the center of mass dynamics is that of a three- 
dimensional anisotropic harmonic oscillator in a rotating 
frame. The characteristic frequencies of this oscillator 
can be calculated from the following classical equations 
of motion determined by the Hamiltonian ( |34| ) (the equa- 
tions of motion and the characteristic frequencies given 
in Ref. |lj are incorrect) 

= -i-R(i) - n x + n x R(t)) .(36) 



One may solve this third order equation for u 2 but the 
formulas are quite complicated. We give here the explicit 
solution only in the simple case, when the angular rota- 
tion vector is directed along one of the principal axes, say 
the z axis, of the potential ellipsoid and the axes x and 
y are directed along the two remaining principal direc- 
tions. In this case the frequency of oscillations along the 
z direction is not modified by the rotation, lo z = Ja^. 
The frequencies of the oscillations in the perpendicular 
directions become equal to 

J2SI 2 + a+± Jal + 8Sl 2 a+ 
= 1 55 • <») 

where a± = a x ± a y . The spectrum of Hcm is discrete 
when the frequencies ui± are real. This is true (cf. @) for 
slow rotations (SI < y/aZ) or for fast rotations (fi > */%), 
where we have assumed for definitness that a x < a y . In 
the intermediate region, the classical oscillations are un- 
stable. This means that for the intermediate values of the 
rotational frequency the spectrum of Hqm is continuous 
and the condensate as a whole falls out of the trap. The 
counterintuitive result that for fast rotations the oscilla- 
tions are stable is due to the action of the Coriolis force. 
The stabilizing effect of the Coriolis force is well known in 
other physical situations. It plays the crucial role in the 
Paul trap and also in the dynamics of nonspreading 
electronic wave packets in Rydberg states, called Trojan 
states |, §. 

The appearance of an unstable region of rotational ve- 
locities between the stable regions of slow and fast rota- 
tions is not limited to the special case discussed above, 
but is a general property of the anisotropic rotating os- 
cillator. This can be seen from the properties of the 
free term in the characteristic equation (p7\). This term 
(taken with the minus sign) is always equal to the prod- 
uct of the squares of the three eigenfrequencies. In the 
absence of rotation this term is equal to In turn, 



The characteristic equation for this system may be ex- 
pressed in terms of six invariant quantities built from A 
and SI and it has the form 



(37) 



for fast rotations this term is again positive since it is 
dominated by SI 2 SI- A- SI. The signature of an unsta- 
ble behavior is the change of sign that indicates that one 
of the characteristic frequencies becomes imaginary. In 
the Appendix we show that for an anisotropic trap this 
change of sign always takes place. Therefore there always 
exists an intermediate region of rotational velocities for 
which the center of mass escapes from the rotating trap. 



V. CONCLUSIONS 

We have made a full use of the fact that for a many- 
body system of interacting atoms in a general harmonic 
potential the dynamics of the center of mass decouples 
from the internal dynamics. In the time-dependent de- 
scription, this leads to the appearance of families of states 
generated by the action of the displacement operator. 
The members of each family are labeled by the solutions 
of the equations of motion for a classical particle in the 
trap. Each family is built on an internal state of the 
system. The internal states are not affected by the dis- 
placement operator. 

For rotating traps — that are often used in BEC ex- 
periments — one may employ the time-independent de- 
scription by using the rotating frame. In that case one 
may study the spectrum characterizing stationary states 
of the system. The decoupling of the center of mass 
Hamiltonian leads to the splitting of the spectrum. The 
spectrum of the full Hamiltonian is a Cartesian prod- 
uct of the center of mass spectrum and the spectrum of 
the internal motion. The spectrum of the center of mass 
Hamiltonian may be discrete or continuous depending on 
the value of the rotational velocity. For sufficiently slow 
rotations the spectrum is discrete because the centrifugal 
force is not strong enough to overcome the trap attrac- 
tion. For sufficiently fast rotations the spectrum is again 
discrete because the Coriolis force stabilizes the center 
of mass motion. In the two stable regions the eigenvalue 



I 

u; 6 - uj 4 (2S1 2 + Tr{i}) + lo 2 ((O 2 ) 2 + 3SI-A-SI - Tr{A} SI 2 + Tr{A} 2 /2 - Tr{A 2 }/2) 
-SI 2 SlA-Sl + Tr{A} SI -A- SI - Sl-A 2 -Sl - Dct{A} = 0. 

I 
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spectrum forms a triplet of ladders with equally spaced 
rungs built on each eigenvalue of the internal Hamilto- 
nian. The spacings are given by the three frequencies, 
solutions of the characteristic equation ([37]). For inter- 
mediate values of the rotational velocity the spectrum is 
continuous that means that the center of mass motion is 
unbounded — the condensate escapes from the trap. 

The spectrum of the center of mass Hamiltonian is not 
influenced by the internal dynamics but the inverse is not 
true. The eigenvalue spectrum of the internal Hamil- 
tonian does depend on the properties of the trap. For 
example, in an exactly soluble model with harmonic mu- 
tual interactions studied previously |(| , the trap potential 
modifies the frequency of the quadrupole oscillations. In 
the present paper we have considered the case of arbi- 
trary two-body mutual interactions but our results are 
also valid in the general case of n-body interactions pro- 
vided the interaction potential is a function of coordinate 
differences only. 

The decoupling of the center of mass motion may 
have direct observational consequences for small conden- 
sates, i.e. for condensates whose physical dimensions are 
small as compared to the extension of the center of mass 
wave function. Such condensates may be produced, for 
example, when interatomic forces are attractive (as in 
lithium). Then, all many-particle correlation functions 
are highly peaked at small distances between the par- 
ticles. In this case, the center of mass wave function 
contains significant (probabilistic) information about the 
position of the whole condensate. In particular, the vor- 
tex lines embedded in this wave function will have exper- 
imentally testable consequences; the probability of find- 
ing the atoms of the condensate close to the vortex line 
is small. Since the center of mass wave function obeys 
the one-particle Schrodinger equation in a trap, the ex- 
plicit solutions of this equation exhibiting various vortex 
structures that were studied in detail in Refs. (l(J and 
lO will become relevant. 



APPENDIX A 

In this Appendix we prove that in a rotating 
anisotropic trap there always exists a range of rotational 
velocities i~i where the center of mass motion is unsta- 
ble. In such a region one of the characteristic frequen- 
cies will be imaginary. As mentioned in the main text, 
this proof is based on the analysis of the characteristic 
equation ((37|). The (minus) free term in this equation 
is equal to the product of the three characteristic fre- 
quencies squared. As seen from (|37j), it is a biquadratic 
polynomial in the absolute value of Jl 

lo\uj 2 lo 2 — £l A n-A-n 

-fl 2 (Tr{A} n- A-n - n-A 2 -nj + Dct{A}, (Al) 

where n is the unit vector in the direction of f2. The 
signature of the transition to an unstable region is the 
change of sign of this polynomial occurring at its zeros. 
The polynomial (AL) has zeros if the discriminant A is 
positive 



(Tr{A} n-A n - n-A 2 n 
ADet{A}n-A-n. 



(A2) 



Without loss of generality, we may assume that a x < 
a y < a z . Then, by rearranging the terms, we can easily 
obtain the form of A that exhibits its positivity 



A = (n 2 x a x (a z - a y ) + n 2 a y {a z - a x ) + n 2 a z (a x - a y )) 
+ An 2 n 2 z a y a z (a z - a x )(a y - a x ) > 0. (A3) 
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The vanishing of A means that there is one double root 
for f2 2 . In that case there is no region of instability. This 
occurs only when the trap is not fully anisotropic but has 
a symmetry axis and rotates around this direction. 



[1] J. J. Garci'a-Ripoll, V. M. Perez-Garcfa, and V. Veksler- 
chik, Phys. Rev. E 64, 056602 (2001). 

[2] C. J. Pethick and L. P. Pitaevskii, Physical Review A 62, 
033609 (2000). 

[3] R. Glauber, Phys. Rev. 131, 2766 (1963). 

[4] M. O. Scully and M. S. Zubairy, Quantum Optics, Cam- 
bridge University Press, Cambridge, 1997. 

[5] I. Bialynicki-Birula, Lett, in Math. Phys. 10, 189 (1985). 

[6] Z. Bialynicka-Birula and I. Bialynicki-Birula, Phys. Rev. 
A 33, 1671 (1986). 

[7] H. Paul, Rev. Mod. Phys. 62, 531 (1990). 



[8] I. Bialynicki-Birula, M. Kalinski, and J. H. Eberly, Phys. 

Rev. Lett. 73, 1777 (1994). 
[9] I. Bialynicki-Birula and Z. Bialynicka-Birula, Phys. Rev. 

Lett. 77, 4298 (1996). 
[10] I. Bialynicki-Birula, T. Mloduchowski, T. Radozycki, and 

C. Sliwa, Acta Phys. Polon. A 100 (Supplement), 29 

(2001). 

[11] I. Bialynicki-Birula, Z. Bialynicka-Birula, and C. Sliwa, 
Phys. Rev. A 61 032110 (2000). 



